Hartree-Fock-Bogoliubov Calculations in Coordinate Space: 
Neutron-Rich Sulfur, Zirconium, Cerium, and Samarium Isotopes 



V.E. Oberacker, A.S. Umar, E. TeranQ and A. Blazkiewicz 

Department of Physics and Astronomy, Vanderbilt University, Nashville, Tennessee 37235, 

(Dated: February 4, 2008) 



USA 



m 
o 
o 

(N 

(D 

ON 



oo 

o 

On 
O 

m 
o 



X 



Using the Hartree-Fock-Bogoliubov (HFB) mean field theory in coordinate space, we investigate 
ground state properties of the sulfur isotopes from the line of stability up to the two-neutron dripline 
^34-52 ^-j particular, we calculate two-neutron separation energies, quadrupole moments, and 
rms-radii for protons and neutrons. Evidence for shape coexistence is found in the very neutron- 
rich sulfur isotopes. We compare our calculations with results from relativistic mean field theory 
and with available experimental data. We also study the properties of neutron-rich zirconium 
^102,104^^^^ cerium (^^^Ce), and samarium C'^^'^^^ Sm) isotopes which exhibit very large prolate 
quadrupole deformations. 



PACS numbers: 21.60.-n,21.60. Jz 



I. INTRODUCTION 



One of the fundamental questions of nuclear structure 
physics is: how many neutrons or protons can we add to a 
given nuclear isotope before it becomes unstable against 
spontaneous nucleon emission? The neutron-rich side of 
the nuclear chart, in particular, exhibits thousands of nu- 
clear isotopes still to be explored with new Radioactive 
Ion Beam facilities 0. Another limit to stability is the 
superheavy element region around Z = 124 — 126 and 
N — 184 which is formed by a delicate balance between 
strong Coulomb repulsion and additional binding due to 
closed shells . Theoretically, one expects profound dif- 
ferences between the known isotopes near stability and 
exotic nuclei at the neutron dripline, e.g. the appearance 
of neutron halos and neutron skins, and large pairing 
correlations. 

There are various theoretical approaches to the nu- 
clear many-body problem. For the lightest nuclei, e.g. 
^^C, an exact diagonalization of the Hamiltonian in a 
shell model basis is feasible 0, Stochastic methods 
like the shell model Monte Carlo approach W, IBI may be 
used for medium- mass nuclei up to A ^ 60. For heav- 
ier nuclei, theorists tend to utilize self-consistent mean 
field theories; both non-relativistic versions B, p,p , ITol 
m m m Ej and relativlstlc versions H M, Ell have 
been developed. As long as the pairing interaction is rel- 
atively weak, it is permissible to treat the mean field and 
the pairing field separately via Hartree-Fock theory with 
added BCS or Lipkin/Nogami pairing. This works well 
near the line of stability Jj. However, as one approaches 
the driplines, pairing correlations increase dramatically 
and it is essential to treat both the mean field and 
the pairing field selfconsistentl'^^ within the Hartree-Fock- 
Bogohubov (HFB) formalism While the HF(B) the- 
ories describe the ground state properties of nuclei, their 
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excited states can be obtained with the (quasiparticle) 
random phase approximation (Q)RPA [l7l llM Il9l| . 

In this paper we study the ground state properties 
of neutron-rich even-even nuclei up to the two-neutron 
dripline. Besides the large pairing correlations already 
mentioned, HFB calculations face another problem in 
this region: not only does one have to consider "well- 
bound" single-particle states (which determine the struc- 
ture near stability), but in addition there are occupied 
"weakly-bound" states with large spatial extent. Fur- 
thermore, because the Fermi energy for neutrons ep ~ 
at the dripline, virtual excitations into the continuum 
states become important for a proper description of the 
HFB ground state. All of these features represent major 
challenges for the numerical solution. 

Traditionally, the HFB equations have been solved by 
expanding the quasiparticle wavefunctions in a harmonic 
oscillator basis |l2l| . This works very well near the line of 
f3— stability because only "well-bound" states need to be 
considered. However, as one approaches the driplines, the 
numerical solution becomes more challenging: in prac- 
tice, it is very difficult to represent continuum states as 
superpositions of bound harmonic oscillator states be- 
cause the former show oscillatory behavior at large dis- 
tances while the latter decay exponentially. On the other 
hand, a direct solution of the HF(B) equations on a 
finite-size coordinate space lattice does not suffer from 
the above-mentioned shortcomings because no region of 
the spatial lattice is favored over any other region: well- 
bound, weakly- bound and (discretized) continuum states 
can be represented with the same accuracy. Therefore, 
the spatial lattice representation has inherent advantages 
for the theoretical description of exotic nuclei. 

Using our recently developed HFB lattice code for de- 
formed nuclei far from stability jl4l |. we have investi- 
gated the ground state properties of the sulfur isotope 
chain, starting at the line of stability {N — 16) up to 
the two-neutron dripline (which turns out to be = 36 
in our HFB calculations). Our calculations show both 
spherical and quadrupole- deformed g.s. deformations; 
in addition, there is evidence for shape coexistence in 
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the very neutron-rich region. In particular, we calculate 
two- neutron separation energies, quadrupole moments, 
and rms-radii for protons and neutrons. Our HFB calcu- 
lations are compared with results from relativistic mean 
field theory and with available experimental data. 

We have also carried out HFB calculations for some 
recently measured heavier systems: among medium and 
heavy nuclei, ^°^7,r {P2 = 0.45(4)) and ^^^Sm (/32 = 
0.46(5)) are among the most deformed isotopes p8j . 
The large deformation could have its origin in the high 
spin down-sloping orbitals near Z = 38,40,62 and N = 
40,64,96. These large prolate deformations at ^"''Zr and 
^^^Sm are confirmed by Hartree-Fock-Bogoliubov calcu- 
lations carried out in the present work. 



II. HFB EQUATIONS IN COORDINATE SPACE 

Recently, we have solved for the first time the HFB 
continuum problem in coordinate space for deformed nu- 
clei in two spatial dimensions without any approxima- 
tions, using Basis-Spline methods The novel fea- 
ture of our HFB code is that it is capable of generating 
high-energy continuum states with an equivalent single- 
particle energy of hundreds of MeV. In fact, early 1-D 
calculations for spherical nuclei |^ and our recent 2-D 
HFB calculations have demonstrated that one needs con- 
tinuum states with an equivalent single-particle energy 
up to 60 MeV to describe the ground state properties 
accurately near the neutron dripline. Moreover, recent 
QRPA calculations by Terasaki et al. suggest that 
one needs to consider continuum states up to 150 MeV 
for the description of collective excited states. It should 
be mentioned that current 3-D HFB codes in coordinate 

space, e.g. Ref. laia, utilize an expansion of the quasi- 

particle wavefunctions in a truncated HF-basis which is 
limited to continuum states up to about 5 MeV of exci- 
tation energy. Alternatively, an expansion in a stretched 
oscillator basis has also been explored p^ . 

A detailed description of our theoretical method has 
been published in ref. _14] ; in the following, we give a 
brief summary. In coordinate space representation, the 
HFB Hamiltonian and the quasiparticle wavefunctions 
depend on the distance vector r, spin projection a = ±i, 
and isospin projection q — ±^ (corresponding to protons 
and neutrons, respectively). In the HFB formalism, there 
are two types of quasiparticle wavefunctions, and 02, 
which are bi-spinors of the form 



(1) 

The quasiparticle wavefunctions determine the normal 
density Pq{r) and the pairing density Pq{r) as follows 
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p.w = - E E -^LM • (3) 

In the wavefunctions, the dependence on the quasiparti- 
cle energy is denoted by the index a for simplicity. 

In the present work, we use Skyrme effective N-N inter- 
actions in the p-h channel, and a delta interactions in the 
p-p channel. For these types of effective interactions, the 
particle mean field Hamiltonian h and the pairing field 
Hamiltonian h are diagonal in isospin space and local in 
position space. 



and 



h{r<jq,r'a'q') = <5,,,. <5(r - r')/i' , (r) (4) 



h{raq,r'a'q') = 5,,,, 5(r - r')/i^,,, (r) . (5) 



and the HFB equations have the following structure in 
spm-space 



{h-^ - A) 



(6) 



with 



(7) 

The quasiparticle energy spectrum is discrete for l-E] < 
—A and continuous for \E\ > —A @. For even-even nuclei 
it is customary to solve the HFB equations for positive 
quasiparticle energies and consider all negative energy 
states as occupied in the HFB ground state. 



III. 



NUMERICAL METHOD 



Using cylindrical coordinates (r, z,(/)), we introduce a 
2-D grid (ra, ^/a) with a — 1, Nr and (3=1, N^- In 
radial direction, the grid spans the region from to Vmax- 
Because we want to be able to treat octupole shapes, we 
do not assume left-right symmetry in z-direction. Con- 
sequently, the grid extends from —z^ax to +Zmax- Typ- 
ically, Z„iax ~ rmax and Ki2- Nr. 

For the lattice representation, the wavefunctions and 
operators are represented in terms of Basis-Splines. B- 
Splines of order M, Bf'[x), are a set {i = 1,...,A/') of 
piecewise continuous polynomial sections of order Af — 1; 
a special case are the well-known finite elements which 
are B-Splines of order M = 2. By using B-Splines of 
seventh or ninth oder, we are able to represent derivative 
operators very accurately on a relatively coarse grid with 
a lattice spacing of about 0.8 fm resulting in a lattice 
Hamiltonian matrix of relatively low dimension. While 
our current 2-D lattices are linear, a major advantage of 
the B-Spline technique is that it can be extended to non- 
linear lattices (e.g. exponentially increasing) which 
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will be particularly useful for problems where one is in- 
terested in the behavior of wavcfunctions at very large 
distances. 

The four components (n = 1,...,4) of the HFB bi- 
spinor wavefunction ^,i(r, z) are expanded in terms of a 
product of B-Splines 

A/', AO 

Mr., = E E Bt'(rc,)Bf'{z0)cii . (8) 



We construct the derivative operators contained in 
the Hamiltonian with the B-Spline Galerkin method 
while local potentials are represented by the collocation 
method [2I The numerical solution of the HFB 

equations results in a set of quasiparticle wavcfunctions 
at the lattice points. The corresponding quasiparticle 
energy spectrum contains both bound and (discretized) 
continuum states. We diagonalize the HFB Hamiltonian 
separately for fixed isospin projection q and angular mo- 
mentum projection fl. Note that the number of quasi- 
particle eigenstatcs is determined by the dimensionality 
of the lattice HFB Hamiltonian. For fixed values of q and 
O, we obtain A - Nr ■ eigenstates, typically up to 1,000 
MeV. 

In ref. 0, |23| we have investigated the numerical con- 
vergence of several observables as a function of lattice 
box size, grid spacing, and maximum angular momen- 
tum projection ^max- In the case of spherical nuclei, 
our calculations have been compared with the 1-D ra- 
dial HFB results of Dobaczewski et al. and indeed 
there is good agreement between the two. Production 
runs of our HFB code are carried out on an IBM-SP 
massively parallel supercomputer using OPENMP/MPI 
message passing. Parallclization is possible for different 
angular momentum states fl and isospins (p/n). 



IV. NUMERICAL RESULTS AND 
COMPARISON WITH EXPERIMENTAL DATA 

In this section we present numerical results of our 
2D-HFB code and compare these to experimental data 
and other theoretical methods. In all of our calcula- 
tions we utilized the Skyrme (SLy4) [l^l effective N-N 
interaction in the p-h and h-p channel, and for the p- 
p and h-h channel we use a delta interaction with the 
same parameter set as in ref. \l4\ : a pairing strength 
of Vq = —170.0 MeVfm^, with an equivalent s.p. en- 
ergy cutoff parameter Smax = 60 MeV. All calculations 
reported in this paper were carried out with B-Spline or- 
der M = 7 and maximum angular momentum projection 

n - ^ 



A. Sulfur isotope chain up to the two-neutron 
dripHne; shape coexistence studies 

The sulfur isotopes [Z ~ 16) have been investigated 
several years ago by Werner et al. "T] using self-consistent 
mean field models: both Skyrme-HF and relativistic 
mean field (RMF) model calculations were carried out 
using a simple heuristic "constant pairing gap" approx- 
imation. Because of the well-known deficiencies of stan- 
dard pairing theory in the exotic neutron-rich region, we 
have decided to re-investigate the sulfur isotope chain, 
starting at the line of stability {N = 16) up to the two- 
neutron dripline (which turns out to be A'^ = 36 in our 
HFB calculations). Based on the above-mentioned ear- 
lier calculations, one may expect a wide range of ground 
state deformations, and in addition there has been some 
evidence for shape coexistence in this region 7]. Be- 
cause our HFB code has been specifically designed 
to describe deformed exotic nuclei, the sulfur isotopes are 
expected to provide a rich testing ground for our calcu- 
lations. 

Within the single-particle shell model, one would ex- 
pect nuclei such as f|S'28 with "magic" neutron number 
= 28 to be spherical. But the mean field theories 
predict, in fact, deformed intrinsic shapes as a result of 
"intruder" states. Furthermore, in some of these isotopes 
shape coexistence has been predicted 't^. All these phe- 
nomena depend strongly on the interplay between the 
mean field and the pairing field with is correctly de- 
scribed the the HFB theory. Furthermore, the neutron- 
rich N !^ 28 nuclei play a crucial role in astrophysics for 
the nucleosynthesis of the heavy Ca-Ti-Cr isotopes 29]. 

In radial (r) direction, our lattice extends from 0—12 
fm, and in symmetry axis (z) direction from —12, -|-12 
fm, with a lattice spacing of about 0.8 fm in the 
central region. Angular momentum projections Q = 
1/2, 3/2, 21/2 were taken into account. 
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FIG. 1: Two-neutron separation energies for sulfur isotopes. 
The dripline is located where the separation energy becomes 
zero. 

Figure n shows the calculated two-neutron separation 
energies for the sulfur isotope chain. The two-neutron 
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separation energy is defined as 

S2n{Z,N) = Eund{Z,N) - Eund{Z,N -2) . (9) 

Note that in using this equation, all binding energies 
must be entered with a positive sign. The position 
of the two-neutron dripline is defined by the condition 
S2niZ,N) = 0, and nuclei with negative two-neutron 
separation energy are unstable against the emission of 
two neutrons. 

The two-neutron separation energies have been calcu- 
lated using various methods: in addition to HFB calcula- 
tions (i.e. selfconsistent mean field with pairing), we have 
also carried out Hartree-Fock calculations with added 
Lipkin/Nogami pairing (HF-I-LN), and we compare our 
results to the relativistic mean field with BCS pairing 
(RMF) calculations by Lalazissis et al. Experi- 
mental data based on measured binding energies '3(J| are 
available up to the isotope *^S. Fig. ^ shows that both 
the HFB and RMF calculations are in good agreement 
with experiment where available but there are dramatic 
differences as we approach the two-neutron dripline: Our 
HFB calculations predict ^"5' to be the last isotope that is 
stable against the emission of two neutrons. By contrast, 
the RMF approach predicts S2n{Z,N) > at least up 
to ^^S. Our HF-I-LN calculations also yield positive 2n- 
separation energies in the mass region investigated here. 
It should be stressed that, on theoretical grounds, neither 
BCS-type nor Lipkin-Nogami type pairing is justified for 
the very neutron-rich isotopes. Furthermore, it is well- 
known that the HF-I-LN method breaks down at magic 
numbers (the pairing gap does not vanish). 
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FIG. 2: Quadrupole moment for protons (in units of barn) 
for even-even sulfur isotopes. 

In Fig. El we show a comparison of the HFB and RMF 
results for the intrinsic electric quadrupole moments of 
the sulfur isotopes. In most cases, the HFB and RMF cal- 
culations show a similar trend: we observe a region with 
predominantly prolate deformation. Note, however, that 
for the most neutron-rich sulfur isotopes, our HFB the- 
ory predicts a prolate ground state whereas RMF theory 
yields an oblate shape. Direct measurements of electric 
quadrupole moments are only available for two of the sul- 
fur isotopes. The data compilation of Stone |^ yields 



intrinsic electric quadrupole moments of Q^^p = — 0.15& 
for and Q'^'^p = +OMb for ^'^S. Our HFB code yields 
Qhfb = +0.1136 and -I-0.203& for these two isotopes. 
The RMF calculations of Lalazissis et al. ^3 give values 
of Qrmf — -1-0.3396 and -f0.1596, respectively. 

In intermediate-energy Coulomb excitation experi- 
ments, the energies and B{E2) values of the lowest ex- 
cited 2+ state were measured for 38,40,42^ jg^ ^^^^ j^j, 
^'^S nil. The analysis of the measured BE{2) val- 
ues in terms of the simple quadrupole-deformed rotor 
model yields the following experimental quadrupole de- 
formations for the even sulfur isotopes ^^^^''S': |/32^^| = 
0.246(16), 0.284(16), 0.300(24), 0.258(36); our cor- 
responding Skyrme-HFB theory results are 02 — 
0.16, 0.26, 0.25, 0.29. Apparently, the HFB resuhs for 
40,42,44^ are in good agreement with experiment, but 
there is a discrepancy in the case of '^^S: Skyrme-HFB 
predicts a less deformed shape than the experimental 
value. There is an even larger discrepancy between exper- 
iment and the RMF calculations which yield an almost 
spherical shape, = 0.054. 

Both HFB and RMF calculations reveal shape coexis- 
tence in this region, with an energy difference between 
the ground state and the shape isomer that is usually 
quite small. For example, in the case of 10 5*32 our HFB 
code yields a ground state binding energy of —362.56 
MeV with a quadrupole deformation of /32 — 0.11, and an 
oblate minimum at /32 = —0.15 which is only 0.49 MeV 
higher than the ground state. The RMF predicts in this 
case a ground state binding energy of —362.97 MeV with 
oblate deformation of P2 = —0.25, and a shape isomer 
with P2 = +0.179 which is located 0.29 MeV above the 
ground state. 

In Fig. |31 we show a comparison of the quadrupole 
moment for neutrons predicted by our HFB calculations 
and the RMF results of ref.0|. In both cases, the general 
trend is very similar to the result obtained for protons. 

2.0 1 — , , , , , , , , , , , — I 
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FIG. 3: Quadrupole moment for neutrons (in units of barn) 
for even-even sulfur isotopes. 

In Fig. 01 we compare the root-mean-square radii of 
protons and neutrons predicted by our HFB calculations 
and the RMF calculations with BCS pairing Near 
the line of /3-stability, the proton and neutron radii are 
almost identical, but as we approach the 2n-dripline, we 



see clearly the development of a "neutron skin" as evi- 
denced by the large difference between the neutron and 
proton rms radii. For example, in 1QS34 our HFB calcula- 
tions yield r„ — 3.935 fm and Vp — 3.364 fm, respectively. 
In general, the RMF calculations predict larger neutron 
rms radii for mass numbers A > 38 than do our HFB 
calculations. 
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FIG. 4: Root mean square radii of protons and neutrons for 
even-even sulfur isotopes. 



B. Strongly deformed neutron-rich zirconium, 
cerium, neodymium, and samarium isotopes 

Recently, triple-gamma coincidence experiments have 
been carried out with Gammasphere at LBNL 28] which 
have determined half-lives and quadrupolc deformations 
of several neutron-rich zirconium, cerium, and samarium 
isotopes. Furthermore, laser spectroscopy measurements 
|32j | for zirconium isotopes have yielded precise rms-radii 
in this region. These medium/heavy mass nuclei are 
among the most neutron-rich isotopes {N/Z « 1.6) for 
which spectroscopic data are available. It is therefore of 
great interest to compare these data with the predictions 
of the selfconsistent HFB mean field theory. 

A comparison of our HFB results and experimental 
data is given in tabled] The theoretical quadrupole de- 
formations of the proton charge distributions agree very 
well with the measured data of ref. . In addition, our 
calculated proton rms-radius for ^°^Zr is in good agree- 
ment with recent laser spectroscopic measurements (see 
Fig. 4 of ref-li^)- Theoretical HFB predictions are also 
given for the neutron density distributions. 

In table ini we give a more detailed comparison of vari- 
ous theoretical calculations for -'^gf S'm. As stated earlier, 
only the HFB theory provides a self-consistent treatment 
of both mean field and pairing properties. In contrast, 
the HF -I- Lipkin/Nogami and the RMF calculations treat 
mean field and pairing as separate entities. By compar- 
ing the HFB and HF-I-LN results in the first two columns 
with the measured binding energies and quadrupole de- 
formations (last column) we find that our HFB lattice 
calculation yields values which are closer to the exper- 
imental data. Table ^ also shows that while the RMF 
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TABLE L Our HFB results for neutron-rich zirconium, 
cerium, neodymium, and samarium isotopes. The first col- 
umn lists the neutron-to-proton ratio N/Z. Subsequent 
columns display quadrupole deformations (32{p), [32{n) and 
rms-radii rp,rn of protons and neutrons. Recent experimen- 
tal data for quadrupole deformations are taken from ref.|28l|. 
The rms-radius for ^'^^Zr was measured in ref.(33|. 

N/Z fe(p) pr^jp) fe(n) rp (fm) r;-^ (fm) r„ (fm) 



Zr 1.55 0.43 0.42(5) 0.43 4.47 4.54 4.65 

Zr 1.60 0.45 0.45(4) 0.45 4.49 4.70 

Ce 1.62 0.32 0.30(3) 0.33 5.01 5.22 

Nd 1.60 0.37 0.36 5.08 5.27 

Sm 1.58 0.38 0.37 5.13 5.31 



theory reproduces the experimental binding energy quite 
well it seriously underpredicts the strong quadrupole de- 
formation measured in ref. In addition, table ^ 
compares theoretical results for other observables such 
as rms-radii for neutrons and protons, Fermi energies 
(A„,Ap), pairing gaps (A„,Ap), and pairing energies 



TABLE II: Ground state properties for efSm. The first two 
columns give the results of our present work in the Hartree- 
Fock Bogoliubov theory (HFB) and in the Hartree-Fock plus 
Lipkin-Nogami pairing theory (HF-t-LN). The third column 



shows RMF theory results 


[I3, and the last column gives a 


comparison with recent experimental dat 


a P5. 






HFB 


HF+LN 


RMF 


exp 


B. E. (MeV) 


-1,290.2 


-1,286.6 - 


1,291.98 


-1,291.9 




0.375 


0.359 


0.292 


0.46(5) 


r„ (fm) 


5.27 




5.376 




Vp (fm) 


5.15 




5.098 




A„ (MeV) 


-5.63 


-5.37 






\p (MeV) 


-9.21 


-8.93 






A„ (MeV) 


0.31 


0.65 






Ap (MeV) 


0.37 


0.75 






Epa^r{n) (MeV) 


-0.96 


-4.90 






Epa^rip) (MeV) 


-1.19 


-4.57 







In Fig. Owe depict contour plots of the density distri- 
butions for neutrons and protons in ^gf^m. The large 
prolate quadrupole deformation is clearly visible. We 
also observe small density enhancements near the cen- 
ter of the nucleus which are caused by the nuclear shell 
structure. 

Fig.Elshows the corresponding pairing density for neu- 
trons and protons; as discussed in ref. 8], this quantity 
describes the probability of correlated nucleon pair for- 
mation with opposite spin projection, and it determines 
the pair transfer formfactor. We can see that most cor- 
related pair formation in ^gfSm takes place outside the 
central region of the nucleus. 
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FIG. 5: Density distribution for neutrons (left) and protons 
(right) in \5|Sm 




FIG. 6; Pairing density distribution of neutrons (left) and 
protons (right) in ^11 Sm 



V. CONCLUSIONS 

In this paper, we have performed Skyrme-HFB calcu- 
lations in coordinate space for several neutron-rich exotic 
nuclei. The coordinate space method has the advantage 
that well-bound, weakly bound and (discretized) contin- 
uum states can be represented with the same numerical 
accuracy. The novel feature of our lattice HFB code is 
that it takes into account high-energy continuum states 
with an equivalent single-particle energy of 60 MeV or 
more. This feature is crucial when one studies nuclei 



near the neutron dripline . 

We have calculated the ground state properties of 
the sulfur isotope chain [Z ~ 16), starting at the line 
of stability {N /Z = 1) up to the two- neutron dripline 
with a ncutron-to- proton ratio of N/Z Ki 2.2. In par- 
ticular, we have calculated two-neutron separation en- 
ergies, quadrupole moments and rms radii for protons 
and neutrons. In comparing our HFB calculations with 
other theoretical methods (RMF with BCS pairing and 
HF-|-Lipkin/Nogami) we find similar results near stabil- 
ity but dramatic differences near the 2n-dripline (see Fig- 
ures 1-4). For example, our HFB calculations predict 
^'^S to be the last isotope that is stable against the emis- 
sion of two neutrons whereas the RMF approach pre- 
dicts S2n{Z, TV) > at least up to ^^S. For "^^S, the last 
even-even isotope for which experimental binding ener- 
gies are available (Fig. 1), the experimental value for 
the two-neutron separation energy is 3.64 MeV, as com- 
pared to our HFB calculation result of 3.49 MeV and the 
RMF result of 5.24 MeV. Both our HFB calculations and 
the RMF calculations of ref.|0 predict the existence of 
shape isomeric states in the neutron-rich sulfur isotopes. 
In Fig. 121 we compare calculated electric quadrupole mo- 
ments. In most cases, the HFB and RMF calculations 
yield prolate quadrupole deformations. However, for the 
most neutron-rich sulfur isotopes, our HFB theory pre- 
dicts a prolate ground state whereas RMF theory yields 
an oblate shape. Shape coexistence is found both in the 
HFB and RMF calculations, with fairly small energy dif- 
ference between the ground state and the shape isomer. 
Specific results are given for f 15*32. A comparison be- 
tween the root- mean-square radii of protons and neutrons 
clearly exhibits the development of a "neutron skin" in 
the neutron rich sulfur isotopes: for example, in 155*34 
our HFB calculations yield neutron and proton rms radii 
of r„ — 3.935 fm and = 3.364 fm, respectively. 

In connection with recent experiments at Gam- 
masphere, we have carried out HFB calculations of 
medium/heavy mass nuclei with N/Z k, 1.6. In par- 
ticular, we have examined the isotopes ^'^^'^^^Zr, ^^^Ce, 
^^^Nd, and ^^^'^^°Sm. The theoretical quadrupole 
charge deformations for zirconium and cerium are in very 
good agreement with the new data. Also, our calculated 
proton rms radius for ^'^^ Zr agrees with recent laser spec- 
troscopic measurements. Table ^ gives a summary of 
these results and presents predictions for two neutron- 
rich isotopes, ^^^Nd and ^^^Sm, for which experimental 
data are expected to become available in the near future. 
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